Experimental bound entanglement in a four-photon state 
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Entanglement [ , ] enables powerful new quantum technologies [ - ], but in real- world imple- 
mentations, entangled states are often subject to decoherence and preparation errors. Entanglement 
distillation [9, 10] can often counteract these effects by converting imperfectly entangled states into 
a smaller number of maximally entangled states. States that are entangled but cannot be distilled 
are called hound entangled [ ] . Bound entanglement is central to many exciting theoretical results 
in quantum information processing [ - ], but has thus far not been experimentally realized. A 
recent claim for experimental bound entanglement is not supported by their data [^^]. Here, we 
consider a family of four-qubit Smolin states [ ] , focusing on a regime where the bound entangle- 
ment is experimentally robust. We encode the state into the polarization of four photons and show 
that our state exhibits both entanglement and undistillability, the two defining properties of bound 
entanglement. We then use our state to implement entanglement unlocking, a key feature of Smolin 
states [ 1. 
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Bound entangled states are important for several rea- 
sons. First, they represent irreversibility in entanglement 
manipulation: they require the consumption of pure en- 
tanglement to be created via local operations and clas- 
sical communication (LOCC), but no pure entanglement 
can be distilled from them via LOCC [11, 17-lP]. Second, 
they constitute a challenge to develop better entangle- 
ment criteria, as there is no standard efficient way to de- 
tect their entanglement [2, 20]. Third, they are central to 
recent breakthroughs regarding channel capacities [ ]. 
Fourth, despite not being distillable, they still constitute 
a resource for quantum teleportation [21], quantum cryp- 
tography [12, 13], and channel discrimination [ ]. Thus, 
bound entanglement is crucial for developing a more com- 
plete picture of the role of entanglement in quantum in- 
formation. 

One of the most elegant and striking examples of 
bound entanglement is the four-party Smolin state [ ], 
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where the subscripts label the parties, and |^^) are 
the two-qubit Bell states. One may understand the 
Smolin state in the following way: A and B share one 
of four possible Bell states, and C and D share the 
same state, but each Bell state is equally likely and un- 
known. The Smolin state is entangled in the sense that 
it does not admit a fully-separable decomposition of the 
form "^ZkPkpA ^Vb^^c^^d^ with pj^ probabilities and 
p^^r]^^r\^^ states. On the other hand, no pure entan- 
glement can be distilled, neither in the form of bipartite 
nor multipartite pure entangled states [ ]. Thus, the 
Smolin state is, by definition, bound entangled. Despite 
its undistillability, the entanglement in the Smolin state 
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is unlockable in the sense that a joint measurement be- 
tween any two parties can enable a pure maximally en- 
tangled state between the other two [ ]. 

It is evident from equation (1) that the Smolin state 
is separable (or unentangled) in the (AB) : (CD) bi- 
partite cut since it can be written in a biseparable form 
p = '^f.PkpAB ^ ^CD' ^^ ^^^ state is symmetric with 
respect to the exchange of any two parties [16, 23], it 
is separable with respect to all three two-two bipartite 
cuts (AB) : (CD), (AC) : (BD), and (AD) : (BC). 
Following the arguments in [ ], one concludes that no 
entanglement can be distilled between any two parties, 
and this excludes also the distillation of three- and four- 
partite entanglement. Similarly, to prove that an experi- 
mentally prepared Smolin state is undistillable, it is suf- 
ficient to show that all eigenvalues remain positive under 
partial transposition (PPT) across all two-two bipartite 
cuts [ ]. 

A recent work reported the production of a pseudo- 
bound entangled state in liquid-state NMR [ ]. They 
demonstrate sufficient control over their system to im- 
plement the transformations that lead to bound entan- 
glement. Yet, they would need to start in a highly pure, 
rather than a highly mixed state, to generate bound, as 
opposed to pseudo-bound, entanglement. Optical sys- 
tems on the other hand can produce highly pure states. 
In Ref [ ] , the authors claimed an optical demonstration 
of a Smolin bound-entangled state. Their data showed 
that their state was entangled. They reconstructed their 
state using quantum state tomography and applied the 
PPT test for distillability and argued that the eigenval- 
ues are all positive or consistent with zero. However, 
their results show five negative eigenvalues over all three 
partial transpositions with a minimum of —0.02 ± 0.02. 
Since the central point in the PPT test is that the min- 
imum eigenvalue has to be non-negative^ their claim of 
demonstrating bound entanglement is not supported by 
their data. 

In any attempt to generate perfect Smolin states, the 
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FIG. 1: Experimental setup to generate four-photon 
Smolin states, a, We generate a family of Smolin states by 
randomly applying unitaries, using two pairs of liquid-crystal 
variable phase retarders (LCRs), to the initial 10^) as 
10^) CD state produced by two down-conversion sources. For 
each party A, B, C and D, the polarization is analyzed 
with a combination of a half- wave plate (HWP), quarter- wave 
plate (QWP) and polarizing beam-splitter (PBS). Birefrin- 
gent crystals (T and T') are used to compensate for the tem- 
poral walk-off inside the down conversion crystals, and bire- 
fringent crystals (c) compensate spatial walk-off in modes 1 
and 3. b, A two-photon interferometer is used to project on 
10) AC for entanglement unlocking between parties B and D. 
The delay At is adjusted for optimum two-photon interfer- 
ence [25]. 



PPT property will be very sensitive to imperfections in 
the state preparation and low counting statistics in the 
experimental data. The main reason is that the partial 
transpose of the density matrix of a perfect Smolin state 
is not full-rank. Introducing a source of white noise leads 
to a full-rank matrix, whose PPT property is more robust 
and thus better suited for experimental demonstration. 
By changing the amount of white noise one generates the 
family of noisy Smolin states 
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where < p < 1 parameterizes the strength of the noise 
applied, and X is the identity matrix. These states are 
bound entangled for < p < 2/3, and fully separable for 
2/3 < p < I [ ]. Entanglement can be ascertained by 
the use of an entanglement witness, that is an observable 
W such that tr(Wr) < for some entangled state r, 
while ti {yVp) > for all separable states p. A suitable 



witness for these states is W = ^— J^i^i ^f ^i^^ i^ i^^]^ 
where cri,cr2 and as are the three Pauli matrices X, Y, 
Z. A short, fully analytic derivation of this witness can 
be found in Appendix A. 

Smolin states can be prepared in the following way. 
Two sources of entangled pairs produce a state of the 
form \^^)ab |(/>+)cD, where |(/)±) = ^ (|00) ± |11)). 
One then applies randomly, but with equal weight, one 
of the rotations a^, p = 0, . . . , 3, with a^ the identity, 
simultaneously to both entangled pairs, i.e., (J^\(I)~^)ab ^ 
cr^|^+)cD- Different levels of white noise can be created 
by choosing a probability, p, where the rotations are ap- 
plied in an uncorrelated fashion. 

In our experiment, we take this approach to generate a 
four-photon Smolin state of the form in equation (2). We 
use spontaneous parametric down conversion to generate 
two photon pairs in four distinct spatial modes, as shown 
in Fig. la. Both pairs are originally in the state |^^), 
encoded in the polarization of the photons, where |0) = 
|H) and |1) = |V) correspond to horizontal and vertical 
polarization, respectively. We implement the unitaries, 
a^, via a pair of liquid-crystal variable phase retarders 
(LCRs) in each of the two sources. For more details see 
the Methods section. 

To characterize our experimentally generated four- 
qubit states, we perform quantum-state tomography on 
the polarized photons with an over-complete set of mea- 
surements. We use an iterative maximum-likelihood al- 
gorithm [2( ] to reconstruct the density matrix. In the 
Methods section and the Appendix B, we provide more 
details on our measurement and state reconstruction pro- 
cedure. In Fig. 2, we show the real and imaginary part 
of the reconstructed density matrices for three differ- 
ent levels of white noise. For p = 0.00, p = 0.49 and 
p = 1.00, the fidelity [ ] with the target state ps{p) is 
(81.52 ± 0.12)%, (96.83 ± 0.05)% and (97.67 ± 0.04)%, 
respectively. The error bars are estimated from Monte- 
Carlo simulations with 500 iterations each. Fidelities for 
the additional noise levels p = 0.25, p = 0.44 andp = 0.75 
are given in Table I. 

In order to test whether our prepared state is bound 
entangled, we directly measure the entanglement wit- 
ness and determine if the reconstructed density matrix is 
PPT. When no white noise is added to our Smolin state, 
only one of these conditions is fulfilled: the state is entan- 
gled but the minimum partial-transpose (PT) eigenvalue 
is negative (see Fig. 3a). The latter is defined as the 
minimum of all eigenvalues of the partially transposed 
density matrix with respect to all two- two bipartite cuts. 
The magnitude of the negative eigenvalue is similar to 
that presented in [ ], which we provide in Fig. 3b for 
comparison. To reach PPT one can add white noise. For 
each level of white noise. Fig. 3a shows the directly mea- 
sured entanglement witness and the minimum PT eigen- 
value of our reconstructed state. As the noise probability 
increases, so do both the witness expectation value and 
the minimum PT eigenvalue. Our state is bound entan- 
gled in the region where the witness is negative while 
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FIG. 2: Experimentally measured density matrices of noisy Smolin states. The upper row represents the real part of 
the measured density matrices for noise levels, p, of (a) 0.00, (b) 0.49 and (c) 1.00 while the lower row contains the corresponding 
imaginary parts. The fidelity with the target state is (a) (81.52±0.12)%, (b) (96.83 ±0.05)% and (c) (97.67 ±0.04)%, and the 
measured witness and smallest PT eigenvalue are, respectively (a) -1.269 ± 0.006 and -0.0273 ± 0.0006, (b) -0.159 ± 0.008 
and 0.0069 ± 0.0008 and (c) 0.985 ± 0.008 and 0.0301 ± 0.0008. From these values we conclude that (a) is entangled but not 
PPT, (c) is PPT but not entangled, but (b) is both entangled and PPT, i.e. bound entangled. 



the minimum PT eigenvalue is non- negative. In partic- 
ular, for p = 0.49 the measured value for the witness 
is -0.159 ± 0.008 (see Table II) and the minimum PT 
eigenvalue is 0.0069 ± 0.0008. These values satisfy both 
conditions required for bound entanglement by a wide 
margin. 

Although the entanglement of a Smolin state (equa- 
tion (1)) is undistillable using only LOCC, it can be un- 
locked between any two of the parties when joint oper- 
ations between the other two parties are allowed [ ]. 
Indeed, by performing a Bell state measurement, two 
parties, e.g. A and 5, can find out which Bell pair they 
share. They can communicate that information to C and 
D^ who will share the same Bell state. In this case, C 
and D will end up sharing a maximally entangled state 
although, as discussed, no entanglement could have been 
distilled between them via four-party LOCC. In the case 
of the family of Smolin states described in equation (2) 
a Bell measurement on any two parties will lead to the 
preparation of a Werner state ["^^J in the other two. The 
resulting Werner state will be entangled for p < |, i.e. 
as long as the Smolin state is entangled. 

To demonstrate entanglement unlocking, we keep the 
noise level at p = 0.49, immediately after obtaining the 
results above, and we feed the photons of parties A and 
C, i.e., one from each source, into the Bell measurement 



depicted in Fig. lb. This performs a projection measure- 
ment of modes A and C onto the Bell state |^~) . An over- 
complete set of two-qubit tomography measurements is 
performed on the qubits of parties B and D^ yielding 
the counts reported in Table III. The reconstructed den- 
sity matrix (see Fig. 4) has a fidelity of (99.45 ± 0.05)% 
with the one expected from the experimental bound en- 
tangled state, assuming a perfect Bell measurement. Its 
tangle is positive, 0.00105 ± 0.00046, and the minimum 
eigenvalue of the partially transposed state is negative, 
—0.0160 ±0.0035; this confirms that we have successfully 
unlocked entanglement between parties B and D. 

We demonstrated experimental bound entanglement, 
for the first time convincingly satisfying its two defining 
criteria: entanglement and undistillability To achieve 
the latter property, we added sufficient white noise to 
clearly fulfill the PPT criterion while still maintaining 
non- separability. Without additional noise, our Smolin 
state is non-PPT by an amount that is relatively small, 
but statistically significant, and we must add a sub- 
stantial amount of noise — almost 50% — to turn it PPT. 
Clearly, it is difficult to achieve the critical PPT condi- 
tion without completely losing entanglement in experi- 
mentally produced Smolin states. Once we achieved the 
preparation of bound entanglement, we demonstrated en- 
tanglement unlocking, realizing all of the conceptually 
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FIG. 3: Experimental tests for bound entanglement. 

Panel (a) shows the measured expectation value for the en- 
tanglement witness (blue diamonds, right axis) and the mini- 
mum PT eigenvalue (red circles, left axis) for various levels of 
white noise. The lines are best fits to the data and the error 
bars correspond to one standard deviation as determined by 
Monte-Carlo simulations based on Poisson distributions for 
the counts. Each eigenvalue shown is the smallest among all 
those calculated from the set of bipartite cuts (AB) : (CD), 
(AC) : (BD), and (AD) : (BC). Using our experimental 
data, we find that our family of generated Smolin states is 
entangled within the shaded region and PPT in the hatched 
region. In the overlapping region, the states are both en- 
tangled and PPT, i.e. bound entangled. In particular, we 
successfully produce a bound entangled state when the noise 
level is p = 0.49; the entanglement witness is —0.159 it 0.008 
and the minimum PT eigenvalue is 0.0069 ±0.0008. For com- 
parison, (b) shows the results from Ref. [ ] , where no addi- 
tional noise was applied. Our state without additional noise is 
entangled but definitely not PPT, with a negative minimum 
PT eigenvalue similar to Ref. [ ]. Our experimental results 
show that a substantial amount of noise is required to turn 
this value firmly positive. 



important characteristics of a Smolin state. Our results 
open the door to applications of bound entanglement in 
experimental quantum information science. 



Acknowledgments 

We acknowledge financial support from NSERC, CFI, 
Ontario MRI ERA and OCE. We thank C. Ryan, O. 
Giihne, R. Laflamme and R. Prevedel for valuable dis- 
cussions. 



Methods 

Two polarized photon pairs are generated using /3- 
BaB204 (BBO) nonlinear crystals phase-matched for 
type-I parametric down-conversion. Each photon source 
consists of two orthogonally oriented, 1mm thick crys- 
tals [29, 30]. The pump laser beam is a 830mW, 
(FWHM= 1.6 nm, centered at 395nm) generated in a 
2mm thick Bismuth-Borate (BiBO) crystal by second 
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FIG. 4: Unlocking of entanglement from a bound en- 
tangled state. By performing a joint measurement on the 
qubits of parties A and C, specifically by projecting them 
on the state |(/)~)ac using the Bell measurement shown in 
Fig. lb, entanglement is unlocked between parties B and 
D. The density matrix of the resulting state shared by B 
and D is reconstructed via an over-complete set of tomogra- 
phy measurements. Here, we use a noise level of p = 0.49 
so that the four-qubit state prepared is bound entangled. 
The reconstructed density matrix of the unlocked entangled 
state (real part on the left, imaginary part on the right) is 
clearly entangled with a negative minimum PT eigenvalue 
of -0.0160 ± 0.0035 and a tangle of 0.00105 ± 0.00046. We 
achieve a fidelity of (99.45 ± 0.05)% with the state expected, 
given the reconstructed four-qubit density matrix of Fig. 2b. 



harmonic generation, from a 2.85W Tiisapphire pulsed 
laser (FWHM= 9.8nm, centered at 790.4nm). Birefrin- 
gent crystals (1mm thick BiBsOe cut at 6> = 152.6° and 
= 0°) are added in modes 1 and 3 shown in Fig. la to 
compensate for the transverse walk-off occurring inside 
the down-conversion crystals. The longitudinal walk-off 
is compensated using a 0.5mm and 2mm quartz plus a 
1mm Qf-BBO crystal before the first SPDC source, and 
a 2mm a-BBO and a 1mm quartz crystal between the 
two sources. The photons are coupled into single mode 
fibers after going through an interference filter with a 
3nm bandwidth centered at 790nm. We initially adjust 
the phase with a quarter- wave plate (QWP) in modes 2 
and 4 such that each source produces |0^). The average 
singles rate produced is 250kHz with an average coinci- 
dence rate of 36kHz, when each analyzer is set to |H) and 
monitoring both outputs of the polarising beam-splitter 
(PBS). 

We use one pair of LCRs per source as depicted in 
Fig. la. Each pair is composed of one LCR that can 
be "off", applying the identity, or "on", acting as an X 
and one LCR acting as either the identity or Z when 
turned off or on, respectively. When they are both on, 
they perform an XZ. For each of the sources and for each 
combination of LCR states. Table IV lists the fidelity and 
the tangle of the state resulting from the LCRs acting on 
the initial |0+) state. At any given time, the state of 
the LCRs is set by a computer using a pseudo-random 
number generator operating at a rate of 10 Hz. 

The polarization of each qubit is analyzed via a half- 
wave plate, and a QWP followed by a PBS (see Fig. 1), 
where we monitor both output modes of the PBS. For 



quantum state tomography, we perform an over-complete 
set of standard polarization measurements |H)/|V) = 
|0>/|1>, |+>/|-> = ^(|H) ± |V)), |R) = ^(|H) + i|V)), 
and |L) = 4-(|H) — z|V)). The counts for each measure- 
ment setting are four-fold coincidence detection events. 
We reconstruct our density matrices using an iterative 
maximum-likelihood algorithm [ ]. A detailed descrip- 
tion of our measurement technique is given in the Ap- 
pendix B. 

In our entanglement-unlocking experiment, we use 
single- mode fibers to feed the photons of parties A and C 
into a two-photon interferometer as depicted in Fig. lb. 
Polarization controllers assure that polarizations |H) and 
|V) are maintained within the fibers. After passing 
through the first PBS in the interferometer, both photons 
are projected on the state |+) using HWPs and PBSs. If 
the time delay Ar is set such that two-photon interfer- 
ence occurs [ ], a coincidence detection event between 
detectors Dl and D2 can only occur if the photons are 
in a state of the form ^ (|HH) + e^^|VV)). By using 
a QWP, we set the phase, x^ such that a coincidence de- 
tection event signals a projection of the two-photon state 
on 1^-). 



check that max^^^ tr(p5'T) < 1/8. We have 



maxtr(p5'r) 
= max tr(p5|a)(a|A |^)((/>|bcd) 

\0i)A\(f>)BCD 

= niax tT{p^s^{\a){a\)^^\(l)){(l)\BCD) 

\a)A\(p)BCD 



1 — 4 min tT{ps\a){a\A 

\a)A\(f>)BCD 



BCD J 



(Al) 



where p^^ denotes the partial transpose of ps with re- 
spect to A. The first equality comes from the sym- 
metry of Ps and the convexity of 5, so that it is 
sufficient to consider pure states \a)A\(l^)BCD in the 
maximization; the second equality from the identity 
tr(Xy') = tr(X^^y'^^); the third equality from p^^ = 
|(cr2 cr^^)l{X-4.ps) {(72 cr^^) and the fact that 
(J2{\c^) {o:\)^(j2 is also a pure state. Finally, the minimum 
in the third line is easily seen to vanish. 

Having established that ps{2/3) is the nearest state 
to ps in 5, according to Theorem 6.1 in [ ] one can 
construct a witness for p^ as W = cqX + psi^/S) — ps, 
with Co = tr(p5(2/3)(/)5-p5(2/3)) = 1/24, so that W ex 
I - Yfi=i crf^ = W. As PS {2/ 3) has fuh rank. Theorem 
6.1 in [31] also assures that W is optimal, i.e., there is 
no witness W^ that detects all the states detected by W 
and some more. 



Appendix A: An entanglement witness for the 
Smolin state 



Appendix B: Quantum state tomography of Smolin 
states 



We consider an entanglement witness W such that a 
negative expectation value (W) = tr(pW) < is suffi- 
cient to exclude that the prepared state p is in the set 
S of mixed states that are a convex combination of pure 
states with one party disentangled from the others, e.g. 
\o^)a^ \(t')BCD' In such a convex combination, the disen- 
tangled party may differ from pure state to pure state. Of 
course, 5 is a superset of the set of completely separable 
states. 

The witness reads W = X — Xli^i ^f^- There are a 
number of ways to single out W as an appropriate wit- 
ness; e.g., in [ ] the derivation was based on the sta- 
bilizer formalism and numerical optimization. Here we 
present a derivation that is completely analytical and 
uses the geometric approach of [oi]. 

It will be sufficient to prove that the state ps{2/3)^ 
which is fully separable [ ] , is the closest to ps in the set 
S. Here, closeness is define d with res pect to the Hilbert- 
Schmidt norm ||X||hs = \/tr(XtX). To this aim, using 
Proposition 5.1 in [ ], one derives that it is sufficient to 



We perform an over-complete set of tomographic mea- 
surements for each level of white noise in order to recon- 
struct the density matrices of the states prepared. The 
measurement for each of the four qubits is chosen from a 
set of six states: |H), |V), |+), |— ), |R), and |L), result- 
ing in an overall set of 6^ = 1296 measurements. Each 
qubit measurement is implemented via a half- wave plate 
(HWP), followed by a quarter- wave plate (QWP) and 
a polarising beam splitter (PBS). By monitoring both 
PBS outputs simultaneously one can reduce the number 
of measurements to 3^ = 81 but this approach requires 
the two outputs of each PBS to have equal coupling ef- 
ficiencies for all analyzer waveplate settings. Even if one 
tried to balance the coupling efficiencies, imperfections 
would lead to systematic errors in the measured expecta- 
tion values, and balancing the coupling efficiencies is very 
susceptible to long-time drifts in the setup. We choose 
a more accurate and less sensitive approach as described 
in the following. Let us assume that we want to measure 
a qubit in the H/V basis. First, we choose the angles 
of the HWP and QWP such that the path transmitted 
through the PBS corresponds to |H), the reflected one 



to |V). We record these counts, and then we repeat the 
measurement but with the waveplates oriented such that 
the transmitted path is |V), and the reflected one is |H). 
By adding the respective counts for these two cases, we 
automaticahy average over any imbalance in the coupling 
efficiencies of the two outputs of the PBS. The approach 
is the same for measuring in any other basis, and for 
each basis we have to average over two settings. For 
four qubits we have to average over 16 settings, resulting 
in a total of 1296 settings, but for each of them we get 
approximately 16 times the number of counts we would 
expect when monitoring only the transmitted paths of 
our PBSs. In other words, using the counts from all PBS 
outputs this way effectively increases the measurement 
time per setting by a factor of 2^, where N is the num- 
ber of two-output qubit analyzers involved. 

Because of the high number of measurement settings 
in each tomographic scan, we take additional precautions 
in order to prevent long-time drifts from influencing our 
measurement results. For each setting we record data for 
5s, and we repeat the full set of measurements 10 times. 
The order in which the settings are measured is random 
for every loop we perform. We add together the counts 
of all loops such that the resulting counts correspond to 
a measurement time of 10 x 5 = 50s per setting. Taking 
into account the additional factor of 16 due to our use of 
all PBS outputs, the effective overall measurement time 
is 800s per setting. 

The method we use for the reconstruction of the den- 
sity matrix is an iterative maximum likelihood technique 
introduced by Jezek et al. [ ]. We compared its results 
with those of another popular technique that poses the 
search for the maximum-likelihood density matrix as a 
semidefinite program [ ], which we solve using CSDP 
(a C Library for Semi-Definite Programming). The lat- 
ter method yields a global optimum for the density ma- 
trix but is in general significantly slower than the iter- 
ative maximum likelihood technique. For comparison 
we performed tomography of the counts measured for a 
noise level of 49% using both techniques. The fidelity 
between the reconstructed density matrices using the 
semidefinite maximum-likelihood technique and the it- 
erative maximum-likelihood algorithm is 1.0000. In case 
of the iterative technique we performed 500 iterations, in 
case of the semidefinite technique we performed 90. Us- 
ing the semidefinite technique we calculate the minimum 
eigenvalue of the partial transpose of the reconstructed 
density matrix with respect to all two-two partitions 
(short: minimum PT eigenvalue) to be 0.007 ± 0.001, 
compared to the 0.0069 ± 0.0008 we get using the itera- 
tive algorithm. These values are in good agreement. 



Appendix C: Minimum eigenvalues and the 
entanglement witness 

Before each loop in the tomographic measurements we 
performed a direct measurement of the entanglement wit- 



ness. Again, we sum the counts of all 10 loops and use 
these results for the witness at each white-noise level. 
To measure the witness we perform 48 measurements. 
These correspond to the 3 x 16 measurements for mea- 
suring the expectation values of af^, for af ^, and for 
af^. The expectation value of the witness is then given 

We give the counts for the witness measurement for 
a noise level of 49% as well as the expectation val- 
ues for the three products of Pauli operators in Ta- 
ble I. Using these results we get an expectation value of 
(W) = -0.159 ±0.008. Table II shows the witness values 
measured for various noise levels, and, for comparison, 
we also provide the witness values calculated from the 
reconstructed density matrices. 

In Fig. 5 we show histograms of all witness values and 
minimum PT eigenvalues for the 500 density matrices 
calculated in the course of the Monte-Carlo simulation 
performed on the counts recorded for a noise level of 49%. 
These histograms further illustrate that the state fulfills 
the criteria for bound entanglement. 

For our entanglement-unlocking results, we provide the 
count rates measured for the tomography of the two- 
qubit density matrix shown in Fig. 4 of the main text in 
Table III. The histograms in Fig. 6 show that the tangle 
is greater than zero, and that the minimum eigenvalue of 
the partially transposed density matrix is negative. The 
histograms comprise the values calculated in 500 itera- 
tions of a Monte-Carlo simulation using the measured 
data. They provide evidence for entanglement between 
the remaining photons, i.e., we have successfully unlocked 
entanglement from an initially bound entangled state. 



Appendix D: Characterization of the 
doAvn-conversion sources 



Depending on the state of the LCRs, different unitaries 
are applied to the states generated by our sources. Both 
sources are aligned to produce \(j)~^) = 4=(|HH) + | VV)) 
states. One of the two LCRs in each source performs 
the identity (ctq) when in the off state, and it performs a 
Pauli X when in the on state. The second LCR performs 
(Jo when off and a Pauli Z when on. The combination of 
the two LCRs determines the overall unitary applied to 
the entangled state. 

In order to test the performance of the two sources, we 
perform two-qubit tomographic scans on the photon pairs 
emitted by the two sources. We do these tomography 
measurements for all four states of the LCRs (off-off, off- 
on, on-off, on-on) for each of the sources. For each set of 
these tomographic measurements we perform a Monte- 
Carlo simulation with 500 iterations, using the measured 
counts as the means of Poissonian distributions; these 
two-qubit measurements were taken immediately before 
the tomography of the Smolin state with 49% noise, and 
are given in Table IV. The values for fidelity and tangle 
are typical for the performance of our down-conversion 
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FIG. 5: Histograms for the entanglement witness and the minimum PT eigenvalue. We used 500 Monte-Carlo 
iterations on the tomography data for p = 0.49 to generate histograms for, a, the expectation value of the entanglement witness 
(W) and, b, the minimum PT eigenvalue. These histograms show that the expectation value of the entanglement witness is 
strictly negative, and that the minimum PT eigenvalue is strictly positive. 
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FIG. 6: Histograms for the tangle and minimum PT eigenvalue of unlocked entanglement. Using the two-qubit 
tomography data for the qubits of parties B and D after projecting the two other qubits on |(/)~), we perform 500 Monte-Carlo 
iterations. Panels a and b show histograms for the tangle and the minimum eigenvalue of the partially-transposed reconstructed 
density matrix, respectively. Here, we use minimum PT eigenvalue in the same meaning as we do for the four-qubit Smolin state 
although, of course, for two qubits there is only one partition with respect to which the partial transpose can be performed. 
Both histograms confirm that entanglement has been unlocked. 



sources. 
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(af *) = 


= 0.4005 =b 0.0043 


(-D - 


= 0.3621 ± 0.0043 



TABLE I: Counts for the measurement of the entanglement witness. The counts are four-fold coincidences measured 
to determine the expectation value of the entanglement witness, (W) = 1 — X^^^i {(yf ), for a noise level oi p — 0.49. We 
integrated over 60s (6s in each of the ten loops) per measurement setting. Because we use all PBS outputs, the effective 
measurement time is 16 x 60 = 960s per setting. 



noise level (in %) 







25 



44 



49 



75 



100 



fidelity 81.52 ±0.12 95.11 ±0.07 96.09 ± 0.06 96.83 ± 0.05 97.40 ± 0.03 97.67 ± 0.04 

min. EV -0.0273 ± 0.0006 -0.0077 ± 0.0008 0.0015 ± 0.0008 0.0069 ± 0.0008 0.0206 ± 0.0008 0.0301 ± 0.0008 

(West) -1.261 ±0.005 -0.673 ±0.007 -0.239 ± 0.008 -0.164 ±0.008 0.419 ± 0.008 1.017 ±0.009 

(Wsum) -1.269 ±0.006 -0.682 ±0.007 -0.253 ± 0.007 -0.159 ±0.008 0.432 ± 0.007 0.985 ± 0.008 



TABLE IL Characterization of the experimentally prepared noisy Smolin states. For various levels of white noise, 
we show the fidelity of the reconstructed density matrix with the theoretically expected density matrix of a noisy Smolin state. 
With "min. EV" we denote the minimum PT eigenvalue. At every noise level, we give two values for the expectation value of 
the entanglement witness. The value we estimate from the reconstructed density matrix is (West), whereas we denote the value 
that we directly measured experimentally (see Table I) as (Wsum)- These values are in good agreement. In the latter case we 
calculate the error using error propagation and assuming Poissonian errors for the measured counts. The values for the other 
quantities (fidelity, minimum PT eigenvalue and (y\^est)) are calculated using the reconstructed density matrix, and the errors 
are estimated via Monte-Carlo simulations with 500 iterations each. 
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3176 


2204 
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TABLE III: Counts for the tAvo-qubit tomography measurement for our demonstration of entanglement unlock- 
ing. We project the qubits of parties A and C on |(/)~) and perform an over-complete set of tomography measurements on 
the qubits of parties B and D. Each row rorresponds to a fixed setting for B, and each column to a fixed setting of D. All 
counts are four-fold coincidence events. We repeated the set of 36 measurements 20 times, and for each measurement setting 
we integrated over 60s. Taking into account the additional factor of 4 due to our use of all PBS outcomes, this results in an 
overall measurement time of 4 x 20 x 60 = 4800s per setting. 





LCR2 


ideal state 


fidelity 


tangle 


LCRl 


source 1 source 2 


source 1 source 2 


cro 


cro 


\4>+) 


95.77 ±0.02% 96.34 ±0.02% 


0.8709 ± 0.0008 0.8685 ± 0.0007 


cro 


cr3 


\4>-) 


96.38 ± 0.02% 96.30 ± 0.02% 


0.8690 ± 0.0007 0.8698 ± 0.0007 


cri 


cro 


\^+) 


96.58 ± 0.02% 96.49 ± 0.02% 


0.8705 ± 0.0007 0.8688 ± 0.0006 


cri 


(73 


IV'-) 


96.12 ± 0.02% 96.51 ± 0.02% 


0.8668 ± 0.0007 0.8664 ± 0.0006 



TABLE IV: Source characterization. For each of the two sources we give the tangle of the state produced and its fidelity 
with the ideal entangled state we want to prepare (|0=^) = ^(|HH) ± |VV)) and |^=^) = ^(|HV) ± |VH)) are the four 
Bell states). Each value given is calculated from the density matrix reconstructed using an iterative maximum-likelihood 
technique on the results of an over-complete set of tomography measurements. The corresponding errors are calculated using 
a Monte-Carlo simulation with 500 iterations. 
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